Fidelity susceptibility of the quantum Ising model in the transverse field 
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We derive an exact closed- form expression for fidelity susceptibility of the quantum Ising model in 
the transverse field. We also establish an exact one-to-one correspondence between fidelity suscep- 
tibility in the ferromagnetic and paramagnetic phases of this model. Elegant summation formulas 
are obtained as a by-product of these studies. 

PACS numbers: 



Quantum phase transitions happen when a small varia- 
tion in an external field can fundamentally change ground 
state properties of a quantum system [l|, [3] ■ They pro- 
vide some of the most striking examples of the richness 
of many-body quantum physics. They can be studied in 
electronic [3| , cold atom [4| , and cold ion [f| systems (see 
Refs. [H, Q for an overview of the field). Traditional con- 
densed matter approaches to quantum phase transitions 
rely on the identification of the order parameter and stud- 
ies of correlation functions A different strategy has 
been recently proposed by the quantum information com- 
munity and is known as the fidelity approach to quantum 
phase transitions [7J, |8( . 

Fidelity is a popular concept in quantum information 
science. It is defined here as the overlap between two 
ground states 

F(g,5) = \(ip(g)\TP(g + 5)}\, 

where \ip{g)) is a ground-state wave-function of some 
Hamiltonian H(g), g is the external field whose varia- 
tion induces a quantum phase transition, while S is a 
small shift of this field. Since the ground states funda- 
mentally change across the critical point, fidelity should 
have a minimum near the critical point [tJ- This sim- 
ple yet powerful observation is the basis of the fidelity 
approach to quantum phase transitions. 

Recent studies suggest that fidelity is an efficient probe 
of quantum criticality In particular, the minimum of 
fidelity near a critical point has been established in sev- 
eral models. The scaling of fidelity with distance from 
the critical point \g — g c \, the field shift 5, and the system 
size N has been shown to encode the critical exponent 
v characterizing the power-law divergence of the correla- 
tion length near the critical point [3, l9Hll| . 

Fidelity has been also shown to play a crucial role in 
quantum phase transitions in quantum fields, where su- 
perpositions of ground states from different phases are 
created [13]. Moreover, fidelity has turned out to be use- 
ful in studies of the dynamics of quantum systems rang- 
ing from simple two-level [l3| to many-body ones (see 



e.g. Ref. [lfj and Sec. V of Ref. [14|]). Its relevance to 
the dynamics of decoherence of a central spin coupled to 
an environment undergoing a q uantum phase transition 
has been described in Ref. [la]. There are two distinct 
regimes where fidelity can be studied. 
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FIG. 1: (color online) Schematic of the phase diagram of the 
Ising model in the transverse magnetic field. 



The first one corresponds to the limit where the field 
shift S is kept constant and the system size N — > oo. In 
this limit one observes Anderson orthogonality catastro- 
phe: Disappearance of the overlap of two ground states in 
the thermodynamic limit [lij . Quite interestingly, while 
the Anderson's seminal paper shows power-law decay of 
the overlap with the system size in a particular model 
that does not undergo a quantum phase transition, there 
is an exponential with the system size decay of the over- 
lap, i.e., fidelity, near the generic quantum critical point 

Qua ( 



(see also Ref. [17]). 



The second regime corresponds to the limit of the field 
shift S taken at the constant system size N. In this 
limit fidelity is close to unity and it can be approximated 
by the lowest-order nontrivial Taylor expansion 

F(g,6)~l- X (g)^, 

where the linear in S term vanishes due to normaliza- 
tion of the ground states. The prefactor in the quadratic 
term, x(flO> is called fidelity susceptibility [l8j]. It has 
been recently intensively studied as a probe of quantum 
criticality [353]. 

We will study fidelity susceptibility in the quantum 
Ising model in the transverse field. This is a paradig- 
matic model of quantum phase transitions: All the basic 
concepts about both equilibrium [l[ and non-equilibrium 
[l9j quantum phase transitions have been tested on this 
model. Its versatile experimental realization shall be pos- 
sible in the nearest future in cold ion setups [20j |. 

The Hamiltonian that we study is given by 
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where g is the magnetic field and N is the number of 
spins. Phase diagram of this model reflects competition 
between the spin interactions and the magnetic field (Fig. 
[I} . The spin interactions try to align spins in the ±x di- 
rection, while the magnetic field polarizes spins along the 
z direction. For large enough magnetic fields the system 
is in the paramagnetic phase, while for small magnetic 
fields it is in the ferromagnetic phase. The two phases 
are separated by the critical point g c = 1. 

To find the ground states of the Hamiltonian ([T]), we 
assume that the number of spins is even and proceed in 
the standard way following notation from Ref. We 
obtain after some calculations 



F(gJ)=Y[cos 



k (g + 6) - 9 k (g) 



. . sin A; 

tan k {g) = -, 

g — cos k 

k = n/N, 3tt/N, . . . , tt - ir/N. 



(2) 



Fidelity susceptibility then reads 

x(g) 



sin k 



4 ^ ( g 2 - 2 5 cosfc + 1) 2: 



(3) 



where the prime denotes the derivative with respect to g. 
Since fidelity susceptibility is symmetric around g = 0, it 
is sufficient to consider magnetic fields g > 0. Further- 
more, it is useful to consider the vanishing magnetic field 
separately. We find with the help of Ref. [2l[ that 



X(0) = 



N 
16' 



From now on, we assume that the magnetic field points 
in the +z direction (g > 0). Moreover, we assume that 
N > 4 excluding the trivial N = 2 case. 

There are at least four options for extracting informa- 
tion out of Eq. ((3]) - note that this sum appears not only 
in the studies of fidelity susceptibility, but also in the sim- 
ilar studies of quantum geometric tensors [22], quantum 
Fischer information [23j ]. and quantum adiabatic evolu- 
tion 24|. The first option is to perform the summation 
numerically, see e.g. Ref. Q- This approach has ob- 
vious limitations. The second option is to expand the 
summand in the Taylor series, see e.g. Refs. 0,113, El- 
This technique is necessarily approximate, it produces 
results whose range of applicability is a priori unknown, 
and it cannot be deployed at arbitrary magnetic fields. 
The third option is to factor out small systems, replace 
y\ by J dk, and calculate the integral, see e.g. Refs. 

[2J, |25| The main drawback of this approach 



is that it produces results whose range of applicability 
can only be guessed. This range of applicability would 
be precisely known, had the "reminders" in the Euler- 
Maclaurin summation formula been studied [26j ] . which 
is complicated and has not been yet done. In particu- 
lar, the replacement of the sum by the integral produces 



meaningless results near the critical point (see the discus- 
sion below). The fourth, ultimate, option is to compute 
the sum exactly analytically, which we will do. 
Our result is based on the following identity 



sin k 



t [sinh 2 (izi/2) + sin 2 (£;/2)]' 



N 2 



cosh 2 (7Vx/2) 



2N 1 



tauh(Nx/2) 
tanh(x) 



(4) 



that we have derived by the end of this paper. Substi- 
tuting x = In g into Eq. Q and comparing the resulting 
expression to Eq. J3|), we find 



x(a) 



N 2 
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I65 2 (g N + 1) I65 2 (9 N + 1) (9 2 - 1) 



(5) 



This result is exact and remarkably simple. In particular, 
it works at and around the critical point, where the most 
interesting physics happens. For example, by rewriting 



-.N 



9 



(1 ' '/ .</ V '■). 



(6) 



we see that Eq. ([5]) is regular at g c 
We notice from Eq. ([5]) that 



9 2 xi.9) 



(7) 



which can be also verified from Eq. ([3]). An analogical 
result was proposed for full fidelity in Ref. [17|. Thus, 
g 2 x(d) is symmetric with respect to the 



ferromagnetic f-> paramagnetic, 
1 

9 -, 

9 



(8) 



mapping. Equation ([7} establishes ferromag- 

netic/paramagnetic duality of fidelity susceptibility: 
All information about fidelity susceptibility is contained 
in one of the phases and can be uniquely mapped to the 
other phase. We will now simplify Eq. ([5]). 

Wc introduce another variable to properly organize the 
following discussion, 



y = Nlng, g = eiq?^j-j 



Its physical meaning is simple, 



V ~ sign (g - 1) 



N 

WY 



where £(<?) is the correlation length of the infinite Ising 
chain in the transverse field g. Note that we used the ex- 
act expression for the correlation length, £(<?) ~ l/ \ lng|, 
instead of the approximate one, £(<?) ~ X/\g — 1|, valid 
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FIG. 2: (color online) Rescaled fidelity susceptibility. Solid 
black line: exact result ©. Red dashed lines: approximations 
near and away from the critical point: Eqs. (|14[) . ([9]), and 
(|12|l without reminders. The inset shows enlargement of the 
central part of the plot. The system size N = 1000. Fidelity 
susceptibility is rescaled by the N(N ~ l)/32 factor. 



near the critical point. The prefactors in both expres- 
sions are of the order of unity and depend on the type 
of spin correlations used to define the correlation length 
(27I . The ferromagnetic/paramagnetic mapping ((5J is 
equivalent to the 

y <-> -y 

mapping. 

Away from the critical point. We define that the 
system is away from the critical point when 

1*1 >1- 

There are two advantages of the reparametrization of x 
in terms of y. First, while \g — 1| is bounded from above 
in the ferromagnetic phase, \y\ can be arbitrarily large in 
both phases making it a good parameter to consider the 
far away limit. Second, as can be seen in Ref. [lH, it is 
also a convenient variable for the actual calculations. 

We start by rewriting fidelity susceptibility on the 
paramagnetic side, g > 1, as 



x(g) 



N 



16 5 2 (g 2 

N 



1) 



l6g 2 (l + g N ) \l + g 



R(g), 

Ng N 



R 




x-R 


< 2e 2 yexp ^ 



(9) 
(10) 
(11) 



Such splitting of fidelity susceptibility into the leading 
term and reminder R{g) is exact. While the leading 
term was known before [1, 0, UK, reminder (fT0|) was 
unknown. Knowledge of the reminder is crucial for de- 
termining where x(ff) can t> e approximated by the leading 



term. Bound (jlll) . valid for y > 1 28], shows that for 
y >• 1 the reminder is negligible (see also Fig. [2]). 

The leading term was known before because it follows 
from the above-mentioned replacement of the sum by the 
integral in Eq. ([3]). Note that without knowing the range 
of applicability of this approximation, one may draw a 
wrong conclusion about singularity of fidelity susceptibil- 
ity near the critical point, i.e., in the limit of g — > g c = 1 
(see e.g. Ref. 22| where one of the diagonal elements of 
the quantum geometric tensor is equal to fidelity suscep- 
tibility that we study). Our result shows why the singu- 
larity is absent: The singularity of the leading term in Eq. 
(|9]) is exactly cancelled by the divergent part of reminder 
(pTO]) . This can be verified using Eq. ([6]). Thus, fidelity 
susceptibility near the critical point is regular rather than 
singular, which we will discuss below. 

Using ferromagnetic/paramagnetic duality ([7]), we 
readily obtain fidelity susceptibility on the ferromagnetic 
side, < g < 1, 

x0-/g) n ,6 s fl(Vfl) 



x(g) 



9" 



16(l-ff 2 ) 



- i?, R = 



sr 



(12) 



and analogical remarks to the ones formulated to dis- 
cuss Eqs. (f9j)-(|llj) apply here. In particular, bound on 

R/(x~ R) is the same as in Eq. (jlip after replacing y 
by — y. Despite differences in the functional form of the 
leading terms on both sides of the critical point - Eqs. 
© and (fl"2")) - the two expressions are in fact two sides 
of the same coin due to symmetry (J7]). 

Near the critical point. Wc define that the system 
is near the critical point when 



\y\ < i. 



(13) 



We start by noting that fidelity susceptibility does not 
have maximum at the critical point g c = 1. This is a 
consequence of ferromagnetic/paramagnetic duality ([7]) 
allowing us to write 



x(g) 



Ha) 
g 2 ' 



hi,,) = /,(- 



From the symmetry of h(g), we see that h(g) has an 
extremum at g c — 1 (as it is not the constant function 
here). This implies that x(g) cannot have maximum at g c 
due to the l/<? 2 factor shifting the maximum to g m < gc- 
To simplify the exact result for fidelity susceptibility 
near the critical point, we Taylor-expand h(g(y)) in y and 
keep the 1/g 2 factor unchanged. Such an approximation 
is welcome for two reasons. First, it is efficient because 
h(g(y)) is an even function of y due to the g ■<-> 1/g 
symmetry of h(g): odd expansion terms vanish. Second, 
it preserves symmetry (0. We get 

7Y(7Y- 1) /, N + 1 {N\ng) 2 
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FIG. 3: Illustration of accuracy of the approximate expres- 
sion for the position of maximum of fidelity susceptibility. 
g m (approx) is given by Eq. (|17[) . while g m (exact) comes from 
numerical studies of the exact expression (JSJ). For small sys- 
tem sizes N the shift of the maximum from the critical point 
predicted by Eq. JT7} is off by a few percent and quickly 
decays with the system size. 



where iV In g = y is a small parameter near the critical 
point (|T3|). We see from bound <jT5j) that for \y\ < 1 
reminder r(g) is negligible; r(g) can be obtained by com- 
paring Eqs. (0 and (|14p. The bound on r(g) is derived 
in Ref. |28| . The lower bound is valid for any y, while the 
upper bound is valid for \y\ < vTO. The closest result to 
Eq. (fT4")l was published in Ref. 23[. Apart from lacking 
the reminder, it suffers from computational errors that 
we list in Ref. El. 



Right at the critical point, Eq. (|T4l predicts 

N 2 N 
X(fc=l)= 32 "32- 



(16) 



This has to be compared to predictions of the scaling the- 
ory of fidelity susceptibility proposing that the leading (in 
the system size) contribution to x(ffc) is proportional to 
./V 2 0, [l(| ■ Note that both the leading and the subleading 
term is exactly captured by Eq. (|TB)) . 

Slightly away from the critical point, we enter the 
regime where the term ~ (iVlng) 2 starts to play a role. 
This term is responsible for shifting the maximum of fi- 
delity susceptibility away from the critical point. Indeed, 
forgetting about the reminder in Eq. (fT4"|) and solving 
equation x'idm) = 0, we obtain (see also Fig. [3]) 



9n 



6 

iV 2 



(17) 



The limit of N ^> 1 was assumed to simplify the result; 
analogical calculation done on the exact expression §5§ 
gives the same result in the considered order. This shows 
that indeed the maximum of fidelity susceptibility is lo- 
cated near the critical point on the ferromagnetic side, 
which we proposed from symmetry ([7]). The maximum 
moves to the critical point when N — > oo. 

Moving further from the critical point, we stay in the 
regime where the reminder is still negligible, while the 
term ~ (iVlng) 2 describes decrease of fidelity suscep- 
tibility from the maximum (Fig. [5]). Equation (|14[) 



shows then how the scaling prediction near the critical 
point, x ~ -^ 2 j gradually breaks down. Finally, for 
|iVlng | = \y\ ~ 1 the reminder becomes non- negligible 
and the crossover to the far away limit begins. 

Derivation of the identity. We start from an iden- 
tity that can be found in Ref. [ill ]. 



E 



sm 2 (k/2) tanh(z/2) 
sinh(;z) 2 



iVtanh 



Nz 



(18) 

valid for summations over the same k as in Eq. © . Mul- 
tiplying both sides of Eq. (|18p by tanh(z/2) and taking 
the derivative of the resulting equation with respect to z, 
we obtain another identity: 



E 

k 



sin 2 0/2) 



t [sinh 2 (z/2) + sin 2 (fc/2)] 
N d 



(19) 



sinh(z) dz 



Nz 

tanh ( — — ) tanh 



(I) 



Next, we multiply both sides of Eq. (JTHJ) by cosh 4 (z/2) 
and again differentiate the resulting equation with re- 
spect to z. We obtain after some additional algebra 



d x ^ 

~dz^ 



sin k 



t [sinh 2 (z/2) + sin 2 (fc/2)] 
4 d 



(20) 



cosh 2 (z/2) dz 



cosh 



Next, we integrate Eq. (|20p over z from to x getting 
Si "' J k y - 2iV 2 + 2iV = 

(21) 



V 

k [sinh 2 (.x/2) 
2iVtanh(iVx/2) 



sin 2 (fc/2)] 
2 cosh(iVx) 



N z 



tanh(x) cosh^(iVx/2)■ 
where the left-hand-side is obtained with the help of 



y 1 

^sin 2 (fc/2) 



Z—¥\J 



iV 2 



(22) 



following from Eq. (|19p. Eq. (|2T|) is equivalent to iden- 
tity Q. 

Summarizing, we have derived an exact closed-form 
expression for fidelity susceptibility of the quantum Ising 
model in the transverse field - the quantity that has been 
intensively "approximately" studied in different contexts 
over the last couple of years @,H[l3, 22-25j|. Moreover, 
we have found an exact symmetry of fidelity susceptibil- 
ity showing that all information about it is contained in 
one of the phases and can be easily mapped to the other 
phase. These results allow for an elegant illustration of 
the fidelity approach to quantum phase transitions in ar- 
guably the most important model undergoing a quantum 
phase transition. They also pave the way for analytical 



5 



improvements of several earlier studies involving similar 
summations, see e.g. Refs. 23|,|2J]. Finally, they should 



stimulate analytical investigations of other exactly solv- 
able models. 
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valid for summations over the same k as in Eq. ((2}. For 
the first sum authors of Ref. [23l ] claim that there is an 
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is incorrect. We got the first sum above from Eq. (|22|) . 
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Supplementary material 



Estimations away from the critical point. We will bound here \R/{\~ R)\- This is conveniently done after 
first replacing g by exp(y/N) 



R 



N- 



exp(2y/iV) - 1 



1 + exp(2y/N) 
X-R " [exp(y/2) + cxp(-?//2)] 2 1 + cxp(y) 

Assuming that y > 1 and N > 2, 



1? 



< A- 



exp(2y/A) - 1 



l + cxp(2y/iV) 
[exp(y/2) + exp(-y/2)] 2 1 + exp(y) 



after applying the standard inequality, \a — b\ < \a\ + \b\. Note that it is sufficient to study only positive y, i.e., to 
focus on the paramagnetic side, thanks to the duality symmetry. 
The subsequent "bounding" proceeds as follows 



N- 



exp(2y/A) - 1 



< A (exp(2y/A) - 1) exp(-y) < (e 2 - l) y exp ( -y + 2 {v Ar N) Q{y - N) 



N 



[exp(y/2) + exp(-y/2)] 

The first step above is obvious. To perform the second step, we use the following inequalities from Ref. [30| 

a r - ¥ < r{a - b)^- 1 for < r < 1, 
a r - b r < r(a - 6)a r ~ 1 for r > 1, 

which are valid for positive and unequal a and b (equality happens for r = 0, r = 1, or a = 6). We substituted a — e 2 , 
6=1, and r = y/N to bound ex.p(2y/N) — 1. 
Then, we bound 



1 + exp(2y/A) 
1 + exp(y) 



< (exp(2y/7V) + l)expH;). 



Next we combine the above results and proceed as follows. 
For 1 < y < N we get 



R 



x-R 



< 2/exp(-y) e - 1 



l + exp(2y/iV)\ 



For y > N we get 



R 



< yexv{-y + 2y/N) 



x-R 

Combining these two bounds we get 

R 



1 1 +exp(-2y/N) 



< 2e 2 yexp(-y). 



< 2yeM-y + 2y/N). 



x-R 



<2e y exp -y 



N 



B(y-N) 



Equality in this bound is reached only at y = oo. 

Estimations near the critical point. We will bound reminder r using the following inequalities 

I : tanh(a;) > 



l + x 2 /3' 
II : tanh(x) < x, 



III : tanh 2 (x) > x 



2 2x 4 



IV : tanh(ar) < x - — + — , 

1 1 x 2x i 

V : rr-r > 



tanh(x) x 3 15 
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All of them are valid for x > of interest in our calculations (equalities hold at x = 0). 

Inequality I follows from the discussion presented in Sec. 3.6.13 of Ref. [3lj . Inequality II is proven by considering 

q(x) = x — tanh(x) q'i x ) — tanh 2 (:c). 

Since q(Q) = and q'(x > 0) > 0, we have q(x > 0) > 0, which establishes inequality II. Inequality III follows from 



, 2x 4 x 6 (9 + 2a; 2 ) 

tanh 2 (x) — x -\ > — - 7T > 0, 

3 3(3 + x 2 ) 2 



where inequality I was employed to bound tanh 2 (a;). Inequality IV is proven by considering 



x 



q(x) = x 1 tanh(x) => q' (x) = tanh (x) — x 2 H — x* > 0, 

3 15 3 

where the bound for q'(x) follows from the proof of inequality III. Since q(0) = and q'(x > 0) > 0, we have 
q(x > 0) > 0, which establishes inequality IV. Inequality V straightforwardly follows from inequality IV. 
To bound reminder 



2N g N 2 g N ~g 2 N(N + l) 



N-l(i +g Nf N-l(l+g N )(gi-l) 6 
we replace g by exp(y/N) and rewrite the resulting expression to the form 

iVtanh 2 (j//2) tanh(y/2) N + ly 2 N 



(\ng) 2 -l, 



2(N-1) (N- l)tanh(y/iV) N 6 2(iV-l) 

Note that the reminder is an even function of y thanks to the duality symmetry. 

We bound the reminder from above by bounding tanh 2 (y/2), tanh(y/2), and tanh(y/N) with inequalities III, IV, 
and I, respectively. It leads to 

y 4 18N 2 - 10 + y 2 N y 4 
r < '■ — < 



720 N(N- 1) ~ N- 140' 

where the last step is valid for \y\ < y/lO. 

Next, we show that the reminder is bounded from below by zero by bounding tanh 2 (y/2), tanh(y/2), and 
l/tanh(y/A r ) with inequalities II, I, and V, respectively. This gives us 



y 4 5N 4 - 207V 2 - 96 / 2 , 

^iloAT^TKT^ 2 )- imN ~ V 2 + 5^- 2 - 6108 ' 



